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FORMING GLOBULAR CLUSTER SYSTEMS SEMI-ANALYTICALLY 
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RESUMEN 

Describimos el escenario para la formacion de sistemas de cumulos globulares en galaxias elipticas utilizando 
un modelo semi-analitico de formacion de galaxias. Trabajando en una cosmologia de tipo ACDM, asumimos 
que los cumulos globulares de baja metalicidad se forman a altos redsliifts en fragmentos pre-galacticos, y que 
la consiguiente fusion rica en gas de estos fragmentos conduce a la formacion de los cumulos globulares de alta 
metalicidad. Comparamos nuestros resultados con datos contemporaneos y estudiamos el caso particular de 
la funcioncs dc distribucion de metalicidad estelar y de los cumulos globulares de la galaxia eliptica cercana 
Ccntaurus A. 

ABSTRACT 

We describe a scheme for the formation of globular cluster systems in early-type galaxies using a semi-analytic 
model of galaxy formation. Operating within a ACDM cosmology, we assume that metal-poor globular clus- 
ters are formed at high-rcdshift in prc-galactic fragments, and that the subsequent gas-rich merging of these 
fragments leads to the formation of metal-rich clusters. We compare our results with contemporary data, and 
look at the particular case of the globular cluster and stellar metallicity distribution function of the nearby 
elliptical Ccntaurus A. 
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1. GENERAL 

Globular clusters (GCs) are long-lived tracers 
star formation, ranging from ancient, metal-poor 
stellar populations (e.g.. Galactic halo GCs) to 
young, solar-metallicity stellar populations created 
during interactions (e.g., Schweizer et al. 1996). Due 
to their relatively "simple" nature (they arc idealised 
simple stellar populations), the ages and mctallici- 
ties of GCs can in principle be unambiguously de- 
termined. Hence, the study of GC systems affords a 
unique insight into the formation histories of galax- 
ies. 

We have investigated the formation of the GC 
systems of elliptical galaxies using the fiducial semi- 
analytic model of Cole et al. (2000). We (Beasley 
et al. 2002a) assume GC formation occurs in two 
modes; in pre-galactic fragments with GC formation 
truncated at high redshift, and during the dissipative 
merging of these fragments. 

With these assumptions, we produced 450 real- 
isations of elliptical galaxy GC systems of over a 
range of halo masses i.e. , \.Ox\Q^^h~^ Mq < Mh < 
1.3 X lO^^h-^MQ. 

Using simple stellar population models (Kurth, 
Fritz-v. Alvensleben & Fricke 1999), we can transfer 
our predictions into the observable plane. We have 
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Fig. 1. The distribution in mean ages for 450 model 
GC systems. Top panel: Distribution in look-back time 
for the metal-rich GCs (dark histogram) and metal-poor 
GCs (light histogram). The metal-poor GC histogram 
has been scaled down by a factor of 6 for display pur- 
poses. Bottom panel: Distribution in formation redshifts 
for the GC sub-populations. 
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Fig. 2. Upper panel: Comparison between the semi- 
analytic model stellar metallicity distribution function 
(solid line), and the observed MDF of old stars in NGC 
5128 from Harris & Harris (2002; shaded histogram). 
Lower panel: Comparison between the model GC MDF 
(soUd line) and the observed NGC 5128 GCs (Held et al. 
2002 private comm; shaded histogram). 

compared our model results to contemporary data 
for the GC systems of elliptical galaxies, and the 
principle results from our study are as follows: 

(i) the total number of GCs (N) scales with host 
galaxy luminosity as N oc L^-^^, rather than the ex- 
pected N (x ° for constant GC formation efficien- 
cies. This scaling is very similar to the observational 
data, and in the model is due to a. A4/L oc L depen- 
dence for ellipticals. 

(a) the formation of metal-rich GCs in low- 
redshift mergers does not significantly alter the Sn 
(number of GCs per unit starlight) of the host ellip- 
tical. This reflects the small gas-fractions of merger 
progenitors at later epochs. 

(Hi) the metal-poor GCs have mean ages of 12 
Gyr, whilst the metal-rich GCs have mean ages of 
9 Gyr (e.g. Fig 1). However, the age range of the 
metal-rich GCs in many individual galaxies is signif- 
icant (5 ~ 12 Gyr) with a small fraction forming at 
the present epoch. 



Our results suggest that metal-rich GCs may di- 
rectly trace the merger history of their host galaxy 
whenever star formation occurs. Therefore, the 
metal-rich GCs are remnants of both star formation 
at later times, and the dynamical assembly of their 
host galaxies, possibly providing a key to probing 
the formation epoch of elliptical galaxies. 

Furthermore, our model provides explicit infor- 
mation on the ages and abundances of both the GCs 
and the "field" stars which aggregate to produce 
the final galaxy. Recently, Harris & Harris (2002) 
have constructed a metallicity distribution function 
(MDF) for the nearest giant elliptical galaxy Centau- 
rus A (NGC 5128) using HST. By selecting model 
ellipticals occupying similar environments to this 
galaxy, we can directly compare its predicted and 
observed stellar content (see Beasley et al. 2002b). 

Such a comparison is shown in Fig 2 for the stellar 
MDF of NGC 5128 and its GC system. The model 
MDFs are qualitatively similar to those observed; 
both model and data have stellar components which 
are predominantly metal-rich ('^ O.SZq), but possess 
a small number of metal-poor stars extending down 
to O.OO2Z0. Early gas infall in the model ameliorates 
the so-called "G-dwarf" problem (e.g. Kauffmann 
1996). 

Interestingly, we find that the model MDFs har- 
bour a greater fraction of stars a,t Z > Zq than the 
observations, yielding a broader (by ~ 0.1 dex), more 
metal-rich MDF. This is possibly a result of the fact 
that these outer-bulge observations are missing some 
of the highest-metallicity stars in NGC 5128. Such 
a comparison facilitates the possible reconstruction 
of the star formation history of this nearby elliptical 
galaxy. 
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